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Abstract 

We present a theoretical proposal for the Herzberg chcuit and controlled accumulation of Berry's 
phase in a chirality-based coded qubit in a triangular triple quantum dot molecule with one electron 
spin each. The qubit is encoded in the two degenerate states of a three spin complex with total spin 
S = 1/2. Using a Hubbard and Heisenberg model the Herzberg circuit encircling the degeneracy 
point is realized by adiabatically tuning the successive on-site energies of quantum dots and tunnel 
couplings across a pair of neighbouring dots. It is explicitly shown that encircling the degeneracy 
point leads to the accumulation of the geometrical Berrys phase. We show that only triangular 
but not linear quantum dot molecule allows for the generation of Berry's phase and we discuss a 
protocol to detect this geometrical phase. 
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I. INTRODUCTION 



As discussed by Herzbergi and Berry^, the wavefunction acquires geometric phased, 
Berry's phase, when it is adiabatically moved along a circuit in parameter space of the 
Hamiltonian, the Herzberg circuit, enclosing a degeneracy point. Since only the topology of 
the circuit determines whether the geometrical phase is accumulated, the Berry's phase is 
insensitive to the effects of interactions between the system and its environment. For this 
reason there is interest in attempting to encode and manipulate quantum information in ge- 
ometric phases, for example holonomic quantum computing^ with generalized non-Abelian 
geometric phased. 

Experimentally, Berry's phase in two level systems has already been demonstrated, in- 
cluding neutron^ and nuclear spins^, superconducting qubitsS and a superconducting charge 
pump^. Preceding the successful experiments with superconducting circuits, the relation 
between Berry's phase of a superconducting circuit and other measurable quantities was in- 
vestigated theoretically^iii, including a theoretical proposal for realizing geometric quantum 
computation with superconducting qubits^^. 

In this work, we demonstrate theoretically the generation of Herzberg circuit and Berry's 
phase in quantum states of a three electron complex in a triangular triple quantum dot 
molecule with one electron spin each within the framework of Hubbard and Heisenberg 
models. The two level system, a qubit, is encoded in the two degenerate states of a three 
spin complex with total spin S = 1/3^^"—. An early proposal to generate geometrical phase 
in degenerate one electron quantum levels of a three-atom system was discussed by Herzberg 
and Longuet-Higgins in Ref. |l in 1963. A triple quantum dot (TQD) molecule studied here 
is related to the three-atom system. We define a two- level system, coded qubit,— by the 
two lowest degenerate levels of a half-filled three electron TQD under an in-plane magnetic 
field. It was shown that the quantum states of coded qubit in a TQD can be manipulated by 
tuning the gate voltages^*^. This opens the possibility described in this work to engineer the 
Hamiltonian to undergo adiabatic and cyclic evolution along the Herzberg circuit resulting 
in accumulation of Berry's phase. Recent experimental"— on the linear TQD have already 
demonstrated the high tunability and coherent manipulation of the many-body quantum 
states in a TQD. However, we show that it is not possible to generate a Herzberg circuit for 
a linear triple quantum dot, only triangular triple quantum dot molecule with control over 
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quantum dot energies and at least one tunneling amplitude is needed. 

The plan of the paper is as follows. In SecJTll we describe the system and the Hamiltonians 
for a triangular and linear TQD. In Sec lIIH we construct Herzberg circuit generating Berry's 
phase in a triangular TQD. We show that it is not possible to construct Herzberg circuit for 
a linear TQD. In Sec JIVt a brief conclusion is given. 



II. THE MODEL 



A lateral TQD is defined by metallic gates on top of a two-dimensional electron gas in the 
{x,y) plane at GaAs/AlGaAs heterojunction with three local minima, capable of confining 
a controlled number of electrons. With one electron in each dot, the extended Hubbard 
Hamiltonian reads, 

3 3 ^3 3 

Hhubb = ^ Eihi„ + X] tij{'B)cl^Cj^ 2 ^ VijPiPj + Uifii^hji + ^ gusSa ■ B, 

(1) 

where pi = hi^ + nj^, Ei is the on-site energy, tij is the tunnel coupling between dot i and 
dot j which acquires the Peierl's phase if field B is applied perpendicular to the closed-loop 
structure of a triangular TQD. Vij is the long range Coulomb interaction between dot i and 
j, Ui is the on-site Coulomb interaction of dot i, g is the g-factor of the host semiconductor, 
/iB is the Bohr magneton, and is the spin of the a-th electron. For the present study, we 
set tij = t = -0.05%*, E, = 0, Ui = U = 2.0Ry*, and V^j = V = 0.5Ry* as the initial state 
of the isolated triangular TQD system. Ry* = 5.97 meV is the effective Rydberg in GaAs. 
For linear TQD case, we set ti3 = and V13 = V/2. In this analysis, we assume Ei and tij 
are independently tunable parameters that will be varied to generate Berry's phase. Fig.[T^] 
and Fig. Pd] show the schematic picture of a triangular and a linear TQD, respectively. 



As discussed in Ref. 
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I, the two arrangements of a TQD lead to two topologically different 
Hamiltonians. The low energy spectrum (4 spin-3/2 and 4 spin-1/2 states) of a half-filled 
TQD is mapped onto the Heisenberg modei^^^^^ with one localized spin in each dot, 

H = ^ JijSi ■Sj + ^ gpsSi ■ B + ^ XijkSi ■ (Sj x S^) , (2) 

i<j i i<j<k 

where the exchange interactions Jj_,'s can be derived from the Hubbard model and expressed 
by microscopic parameters as. 



U-V + {E,~ E,) U-V-iE,-E^)J 
The coefficient Xijk for the chirahty operator in Eq.(|2]) is non-zero only for a triangular 

TQD in the presence of a perpendicular magnetic field and is much smaller than Jij . 

Since the Heisenberg Hamiltonian commutes with the total 5*^^ operator of the system, 

we focus on the Sy = —1/2 subspace. We further assume that an in-plane magnetic field By 

has been applied to separate the Sy = —1/2 and Sy = 1/2 subspaces by the Zeeman energy. 

In the Sy = —1/2 subspace, a resonant triangular TQD with all Jij = Jq has the following 

three eigenstates^'^ii^ , 

I g+) = i= (Ut) + m) + Ul)) , (4) 

I g_) = i= [m) + e-^^ m) + e-^^ Ul)) , (5) 

\s,/2) = ^m)+m)+ni))- (6) 

where spin configurations, like (lit); label the electron spins in quantum dots from spin up 
in dot 1 (rightmost) to spin down in dot 3 (leftmost). The two chiral states, and 1^7-), 
constitute the two levels of a chirality-based coded qubit in a triangular TQD. The coherent 
manipulation of a coded qubit is achieved by tuning the exchange interactions Jj/s. In the 
basis of the {|q'±)}, the Heisenberg Hamiltonian with arbitrary Jjj's reads, 

l(2J2s-Jn-Jl2)-^4iJn-J: 



Tjqb 



'12 J 



K2^23 - ^13 - J12) + ^^(^13 - J12) 



(7) 



This matrix shows that the chirality-based coded qubit manipulated by Jj^'s is equivalent 
to spin-1/2 particle under an effective magnetic field in the x — y plane. 

For the case of a linear TQD where J12 = J23 = <^ and J13 = 0, the three eigenstates of 
the system are^*^*^ 

|Lo) = -^(Ut)-m)), (8) 
|Li) = -^(^it)-2m)+m)), (9) 

i53/2) = ^(ut)+m)+m))- (10) 

For the linear TQD, the two Jacobi states, | Lq) and | Li) constitute the levels of a coded 
qubit in a linear TQD. The two Jacobi states used for the coded qubit may be physically 



distinguished by the joint spin states on dot 1 and dot 3. For | Lq), the spins in dot 1 and 
dot 3 form a spin singlet; while the spins in dot 1 and dot 3 form a linear combination of 
spin triplets with Sy = —1 and S'j^ = in | Li). Similar to the coded qubit in triangular 
TQD, we present the Heisenberg Hamiltonian of a linear TQD with arbitrary Jj^'s in the 
basis of {|Lo(i))}, 

|( Ji2 + J23) ^ (-^23 - -^12) ^^^^ 

(-^23 — J12) ~\ {Jl2 + <^23) 

Different from the chirality-based coded qubit, the Jacobi states-defined coded qubit in a 
linear TQD is equivalent to a spin-1/2 particle under an effective field in the x — z plane. 
Finally, we remark that this Jacobi basis of a linear TQD also diagonalizes resonant trian- 
gular TQD Hamiltonian. In a triangular TQD, the ground state is doubly degenerate, so 
the two Jacobi states, | Lq) and | Li), are related to the chirality states | q±) via a unitary 
transformation. 



TTQb 



III. GENERATION AND DETECTION OF BERRY'S PHASE 

In the previous section, we presented the Hamiltonian for a coded qubit in both chirality 
basis (triangular TQD) and Jacobi basis (linear TQD), and compared the coded qubit to a 
spin-1/2 particle, whose Hamiltonian, if(R), is a function of an effective field, R = (X, F, Z), 
expressed in terms of exchange couplings Jij. By adiabatically varying the direction of the 
magnetic field around a closed circuit that encircles the diabolical pointr^'^'^^in the parameter 
space of R, the system undergoes a cyclic evolution and accumulates Berry's phase. The 
diabolical point is the point at which the two-level system is degenerate in the parameter 
space of R. 

To generate Berry's phase, we need to vary exchange couplings to rotate the effective 
magnetic field R. We note that the real and imaginary part of the off-diagonal matrix ele- 
ment in Eq.(I7j) can be independently set to either positive or negative value, if we assume 
the capability to independently control all three exchange interactions. Thus, we can engi- 
neer the effective field R for a chirality-based coded qubit to adiabatically traverse a closed 
circuit that encircles the origin in the parameter space. The next task is then to tune the 
exchange interactions with experimentally accessible quantities, such as Ei and tij. Since 
we need to independently control three exchange interactions, we need to vary at least three 
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variables. Experimentally, it is usually more desirable to tune the on-site energies than the 
tunnel couplings. So the simplest attempt is to vary just the on-site energy of each dot. 
However, this simple attempt fails, because each exchange interaction Jij depends on the 
on-site energies of dot i and dot j through the energy difference AEij = Ei — Ej. We note 
a constraint AE12 + AE23 + AE^i = exists for these three energy differences. Therefore, 
we have actually two degrees of freedom with the three on-site energies. We need to tune 
at least one tunnel coupling. For instance, the effective field can be rotated in a perfect 
circle with a radius J by tuning Ei, E^, and tis, while E2 = E, t23 = tu = t are being held 
constant. Furthermore, we select values of Ei, E^, and tis at each point on the circle such 
that = J is obeyed at all times, and 7 is a constant. Given these conditions, we need 
to tune three parameters according to 



and Fig.[2]D] present the variations of the AE12, AE23 and ^13 as functions of 9, and Fig. [2}:;] 
shows the variations of Jij in response to the changes in the Hubbard parameters. Fig.[2li] 
shows the Herzberg circuit with radius J in the parametric space R for the effective spin- 1/2 
model. Although the circuit presented in Fig. [2] indeed encloses the diabolical point in the 
parameter space, it is based on the Heisenberg model for an effective two-level system. The 
two lowest states in the Sy = —1/2 subspace constitute a two-level system modelled by 
Heisenberg Hamiltonian only when on-site Coulomb repulsion is the dominant energy scale. 
If we examine the scale of variations for A£'i3 and AE23 in Fig. [2], the variations are as high 
as 0.58 U at certain points on the circuit. The validity of the Heisenberg model becomes 
questionable over the course of transporting the system around the circuit. Thus, we expect 
that a more realistic circuit, which minimizes the variations of AEij in order to ensure the 
validity of Heisenberg model, would require to tune more than just three variables. For 
instance. Fig. [3^] and Fig. [3)d] present another circuit, which vary again on-site energy Ei 
and E3 but all ttj in order to produce identical Jij as shown in Fig. [St] and also the same 




(12) 



(13) 



(14) 



where a 



1/4, l3 = V^/A, and 6 is the accumulated angle on the closed circuit. Fig. [2^] 
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circle in Fig. [2]i]. By moderately varying ti2 and ^23, we observe that the variations of AEij 
are significantly suppressed as shown in Fig.[3K]. Fig.[l^] and Fig.[DD] compares the energy 
gap between the lowest three levels in the Hubbard model and the Heisenberg model in the 
Sy = —1/2 subspace for circuit presented in Fig. [2] and Fig. [5], respectively. As shown, the 
second circuit which varies all three tunnel couplings provides a better agreement between 
the Hubbard and Heisenberg model. For many tasks, such as manipulating the quantum 
state of a coded qubit via dynamical phases, a fast tuning of on-site energies is preferred 
to the tuning of tunnel couplings. However, for the adiabatic accumulation of geometrical 
phases the gating operations can be done at a slower rate, with the tuning of the tunnel 
couplings more favourable in this case. 

The accumulated Berry's phase, 4>±{C), for the coded qubit level | q±), after one round 
along the closed circuit C enclosing the origin is simply (f)±{C) = ^i^Q{C), where Q{C) is 
the solid angle subtended by the closed circuit C with respect to the origin of the parameter 
space. For the chirality-based coded qubit, the effective field R is restricted to lie in the x — y 
plane as implied by Eq.(I7]). The solid angle subtended by any closed circuit in the x — y plane 
of the parameter space is either or tt, and this result depends solely on whether the closed 
circuit encircles the origin (diabolical point) or not. Fig. [5] shows the numerical computation 
of accumulated geometrical phase for the coded qubit along the Herzberg circuit presented 
in Fig. [2]. After one round on the circuit, the state accumulates a phase of vr in agreement 
with the theory. To get Fig. [5], we simulate the time evolution, ip{9), of the coded qubit by 
varying the on-site energies and tunnel coupling according to Eq.( |T2|) - Eq.l fT^ . We use the 
angle 6, which denotes the fraction of the Herzberg circuit, as a time variable. Next, we use 
the definition of geometrical phased, 

^jJ{e) = Te-i/^'''^"(^(^))e^^('') l^(RW)), (15) 

where i?„(R(^)) and | n(R,{6))) are the n-th eigenenergy and eigenfunction for the Hamil- 
tonian with parameter R. We remark that | n(R{9))) should be chosen to be single-valued 
and the phases of the eigenstates at different R should be continuously differentiable with 
respect to R. By Eq.( fT5|) . one may extract the accumulated Berry's phase in a numerical 
calculation. 

Let us now turn to the experimental observation of Berry's phase using quantum inter- 
ference. The restriction of having an effective in-plane magnetic field R = {X, Y) makes 
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the experimental probing of Berry's phase difficult because we have very limited interfer- 
ence between components of the wavefunction. Therefore, it is essential to have a closed 
circuit that lies above the x — y plane so the solid angle subtended by the circuit is given 
hj Q = 27r(l — cosG), where cos0 is the direction cosine of the effective field along the z 
direction. To generate an effective field along the z direction, we need to apply a magnetic 
field perpendicular to the triangular TQD. The perpendicular field turns on the chiral term 
in Eq.(|2]). Since the chiral states are eigenstates of the chirality operator, the chirality term 
acts as a operation in the computational space. The coefficient, Xijk, attached to the 
chirality operator is ~ t^/U'^. If we take t = —0.05 and U = 2.0, the values used to gen- 
erate Fig. [3], then Xijk is about two to three orders of magnitude smaller than Jj/s. It is 
important to realize that Xijk, derived from third order perturbation theory, also depends 
on the TQD parameters such as tij, and differences of on-site energies through terms like 
{u-v+AE- ■){u-v+AEki) • Therefore, as the system is adiabatically transported on a closed cir- 
cuit by varying the gate voltage, the magnitude of this effective z field will oscillate in its 
magnitude. Furthermore, an estimated change of the effective z field over the course of one 
complete circuit is two to three orders of magnitude smaller than due to the Ji/s. 

With all the necessary ingredients in place, we now describe a possible procedure to 
experimentally detect Berry's phase in a triangular TQD. First, we prepare the coded qubit 
in a linear superposition of the form ^ | g+) + I Q'-) by having an effective field R 
along the y direction. Next, we follow Eqs. (fT2|) - (fT4|) to tune Ei, E^, and tis under the 
perpendicular magnetic field to accumulate the geometrical phase over one closed circuit 
lying above the x — y plane. Over the course of transporting the system around the circuit, 
the system acquires both dynamical and geometrical phases. Therefore, we next have to 
perform a spin-echo^i^i^ procedure to eliminate the dynamical phases. For this spin-echo 
procedure, we perform a NOT operation on the coded qubit to flip the two chiral states. 
Then we transport the system around the same circuit in the opposite direction. In this 
way, the dynamical phases will add destructively while the geometrical phases will add 
constructively. At the end of the second round along the circuit, the system now has only 
Berry's phase left. The final task is to perform quantum interferometry to extract Berry's 
phase from the system. At this stage, it is important to recall the Jacobi states, | Lq) 
and I Li), which can be distinguished by the joint spin states on dot 1 and dot 3. We 
remark that experiments^^i^i have already demonstrated that the joint spin states between 
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two neighbouring dots can be measured by the spin blockade phenomenon. If the charge 
detection shows that we can make a transition from (1, 1, 1) to (2, 1, 0), then we have a | Lq) 
state in the triangular TQD. Therefore, by simply projecting the chirality-based coded qubit 
onto the Jacobi states, the probability of detecting | Lq) state is given by 



As mentioned in the previous section, the two bases are related by a unitary transforma- 
tion. The measurements in the Jacobi basis thus implicitly provide the needed quantum 
interference to extract Berry's phase. 

We now focus on the case of the coded qubit in a linear TQD. We note that the diagonal 
element Hn of the Hamiltonian, Eq. fllip . is always positive because all Jij > 0, a condition 
dictated by the microscopic details given in Eq.dH]), in which U is the largest energy scale. 
This implies that when we map the coded qubit in a linear TQD onto an effective spin- 
1/2, the effective magnetic field component in the z direction is always positive. Hence 
effective magnetic field in the z — x plane cannot complete a circuit enclosing the orgin of 
the parameter, i.e., effective magnetic field, space. Hence we conclude that it is not possible 
to generate Herzberg circuit and Berry's phase with a coded qubit in a linear TQD. 

IV. CONCLUSION 

In summary, we presented a theoretical proposal for the Herzberg circuit and controlled 
accumulation of Berry's phase in a qubit encoded in the two degenerate chirality states of a 
three spin complex with total spin S = 1/2 in a triangular triple quantum dot molecule with 
one electron spin each. Using a Hubbard and Heisenberg model the Herzberg circuit encir- 
cling the degeneracy point is realized by adiabatically tuning the successive on-site energies 
of quantum dots and tunnel couplings across pairs of neighbouring dots. It is explicitly 
shown that encircling the degeneracy point leads to the accumulation of the geometrical 
Berrys phase. We show that only triangular but not linear quantum dot molecule allows for 
the generation of Berry's phase and we discuss a protocol to detect this geometrical phase 
in interference experiment relying on spin blockade spectroscopy. 




(16) 
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FIG. 2: (Color online) (a) Variations of AE12 (blue curve) and AE23 (red curve) over one round 
on the closed circuit. Note that the variations is comparable to U, the largest energy scale of the 
Hubbard model, (b) Variation of ti3 over one round on the closed circuit, (c) The values of Jij, 
computed with Eq.(l3]), over one round on the closed circuit. J12 is red, J23 is blue, and J13 is green. 
This circuit is generated under the constraint that J13 is held constant, (d) Parametric plot of the 
closed circuit itself in the parameter space R for the coded qubit (effective two-level system) . The 
X and y components of the circle are related to the Jij via the real and imaginary components of 
the off-diagonal matrix element in Eq.d?]). 
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FIG. 3: (Color online) (a) Variations of AE12 (blue) and AE23 (red) over one round on the closed 
circuit. The variations of AEij are significantly suppressed when compared to Fig. [2^]. This is 
because we vary all tij. (b) Variation of Uj over one round on the closed circuit. ti2 is blue, ^23 is 
red, and tis is green. The values of tij are chosen specifically to reproduce the same Jij in Fig.[2h] 
and reduce the variations of AEi. 
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FIG. 4: (Color online) (a) The energy gap £2 — ei (higher ones in the plot) as well as ei — eo (lower 
ones in the plot) in both Hubbard (dotted curve) and Heisenberg (solid curve) model for the circuit 
presented in Fig. [2]. is the i-th eigenenergy of the system, (b) The same energy gaps in both 
Hubbard (dotted) and Heisenberg (solid) model for the circuit presented in Fig. [3]. In both (a) 
and (b), the green, solid curve which represents the energy gap between the two coded qubit level 
is a constant over the closed circuit. This is because we transport the coded qubit on a constant 
energy surface as implied by the parametric plot of {B^, By} for the coded qubit (effective two-level 
system) in Fig.[2]i]. 




FIG. 5: Accumulation of Berry's phase by the coded qubit level, | At the end of the Herzberg 
circuit, 6 = 2tt, the accumulated phase is vr, in agreement with theory. 
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